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1.0  INTRODUCTION 

The  Defense  Advanced  Research  Projects  Agency  (DARPA)  is 
interested  in  exploiting  balloon-borne  communication  and 
surveillance  capabilities  for  military  applications.  Balloons 
offer  several  unique  capabilities  when  high-altitude  balloon 
technologies  are  integrated  with  payloads  employing  the  power- 
efficient,  light-weight  electronic  technologies  available  today. 
In  order  to  establish  communication  and  surveillance  systems 
design  requirements  and  deployment  schedules,  DARPA  must  analyze 
upper-atmosphere  drift  pattern  simulation  results  for  a  variety 
of  balloon  payload  concepts. 

Coleman  Research  Corporation  (CRC),  in  support  of  DARPA 
through  a  Small  Business  Innovative  Research  (SBIR)  Phase  I 
effort,  has  demonstrated  the  technical  feasibility  of  predicting 
high-altitude  balloon  drift  patterns  using  a  digital  computer 
simulation.  CRC  realized  this  ultimate  goal  by  completing  each 
of  four  contract  objectives;  performing  a  literature  survey, 
developing  and  integrating  atmosphere  and  balloon  models  into  a 
Balloon  Drift  Pattern  Simulation  (bdpS),  developing  drift  pattern 
output  formats,  and  exercising  the  BDPS  to  produce  and  analyze 
balloon  drift  patterns.  Drift  patterns  were  produced  for  two 
missions  which  were  specified  by  DARPA.  The  first  mission 
employs  a  zero-pressure  balloon  which  floats  at  an  altitude  of 
70,000  feet  over  24  hours.  The  second  mission  uses  a  super¬ 
pressure  balloon  which  floats  at  120,000  feet  for  a  one-year 
period . 

The  schedule  which  CRC  followed  during  this  six-month  effort 
is  shown  in  Figure  1.0-1.  This  figure  shows  the  sequence  of 
performance  of  the  various  tasks.  The  overlapping  schedule  of 
some  of  the  tasks  indicates  that  mutual  dependencies  of  those 
tasks  required  that  those  activities  be  performed  in  parallel. 

CRC  communicated  throughout  the  contract  performance  period 
with  DARPA  personnel  and  with  DARPA-appointed  contacts.  Through 
this  communication,  several  technical  guidelines  were  received. 
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The  primary  sources  of  information  were  Lt.  Cdr.  Larry  Epley  and 
Lt.  Cdr.  Ted  Krai,  both  of  the  Naval  Technology  Office  of  DARPA 
in  Arlington,  Virginia,  and  Les  Reading,  an  employee  of 
Scientific  Research  Laboratory  in  Santa  Maria,  California.  Lt. 
Cdr.  Epley  described  the  U.S.  Navy's  mission  for  using  the 
balloons  and  requested  that  CRC  model  balloon  ascent  and  descent 
in  addition  to  the  constant-altitude  drift  modeling  originally 
proposed  by  CRC.  Lt.  Cdr.  Epley  also  requested  that  CRC 
concentrate  on  the  24-hour,  zero-pressure-balloon  mission  at  the 
expense  of  the  1-year,  super-pressure-balloon  mission.  Lt.  Cdr. 
Epley  and  Les  Reading  jointly  approved  output  formats.  Lt.  Cdr. 
Krai  specified  six  launch  sites  to  be  used  in  the  production  of 
24-hour  drift  patterns  and  one  launch  site  for  the  production  of 
a  1-year  drift  pattern.  Finally,  Les  Reading  provided  CRC  with 
some  zero-pressure  balloon  design  and  performance  information. 
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2.0  TECHNICAL  APPROACH 


2.1  LITERATURE  SURVEY 

CRC  used  several  resources  in  surveying  the  available 
literature  on  balloon  and  atmospheric  modeling.  In  addition  to 
personnel  at  DARPA,  CRC  initiated  contacts  with  engineers  at 
NASA's  Balloon  Projects  Branch  at  NASA  Goddard's  Wallops  Island 
Flight  Facility,  researchers  at  Texas  A  &  M  University,  and  staff 
at  the  National  Scientific  Balloon  Facility  (NSBF)  in  Palestine, 
Texas.  These  contacts  provided  information  which  guided  CRC  to 
specific  research  reports  which  had  some  relevance  to  BDPS  design 
and  construction. 

Having  identified  specific  reports  for  collection,  CRC 
utilized  the  Redstone  Scientific  Information  Center  (RSIC)  in 
Huntsville,  Alabama.  Many  of  the  pertinent  documents  were  found 
to  be  readily  available.  Other  documents  were  requested  either 
through  RSIC  or  through  the  previously  mentioned  contacts.  As  is 
common  with  surveys  of  this  nature,  the  specific  documents  which 
CRC  set  out  to  find  referenced  many  other  documents,  some  of 
which  also  proved  to  be  useful.  In  addition  to  chaining  through 
technical  journals,  CRC  also  requested  RSIC  to  perform  keyword 
searches  of  on-line  literature  databases. 

A  number  of  scientific  ballooning  reports  provided  CRC  with 
background  information  for  BDPS  design,  construction,  and 
operation.  The  single  most  useful  volume  which  served  this 
purpose  was  the  Scientific  Ballooning  Handbook  (ref.  IJ.  This 
work  contains  discussions  of  the  history  of  ballooning,  the 
physics  of  balloon  flight,  balloon  design,  lifting  gas  selection, 
and  operational  considerations. 

Kreith  (ref.  2]  provides  a  detailed  perspective  on  the 
thermal  design  factors  which  affect  a  scientific  balloon's  design 
and  performance.  He  presents  an  analysis  of  the  energy  balance 
for  the  balloon  system  by  looking  first  at  the  kinematics  of  the 
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balloon  system,  then  at  the  thermodynamic  interaction  between  the 
balloon  and  the  atmosphere,  and  finally  at  the  interaction  of  the 
balloon  fabric  a'.d  the  lifting  gas. 

Carlson  and  Horn  developed  a  balloon  model  which 
incorporates  the  dynamic  and  thermodynamic  aspects  of  high- 
altitude  balloon  flight.  Their  work  is  presented  in  two  volumes, 
the  first  of  which  describes  their  analytical  approach  {ref  3]. 
The  second  document  (ref.  4)  provides  the  details  of 
implementation  of  the  model,  including  FORTRAN  source  code  and 
sample  output  from  execution  of  the  program.  The  source  code 
contained  in  the  latter  of  the  two  documents  proved  to  be  very 
helpful  during  the  development  of  BDPS  because  it  provided 
pertinent  modeling  examples.  It  should  be  noted  that  the  NSBF 
and  NASA's  Balloon  Projects  Branch  both  report  that  they 
currently  use  this  program  (with  minor  modifications)  for  their 
in-house  simulation  purposes. 

A  particularly  informative  report  by  Dwyer  (ref.  5)  was 
discovered  through  one  of  the  on-line  searches  performed  by  RSIC. 
Dwyer  categorizes  and  critiques  existing  aerodynamic-thermo¬ 
dynamic  models  for  predicting  the  vertical  motion  of  zero 
pressure  balloons.  The  report  identifies  deficiencies  of  some  of 
the  existing  modeling  approaches  and  proposes  corrections  and/or 
improvements.  This  report  is  an  excellent,  in-depth,  technical 
review  of  balloon  modeling  practices.  However,  the  factors 
discussed  in  the  report  far  exceed  the  level  of  detail  required 
in  the  present  effort  to  develop  a  balloon  drift  pattern 
simulation . 

The  flight  simulation  models  and  programs  surveyed  model  the 
vertical  motion  of  scientific  balloons  yet  none  were  found  which 
also  address  the  lateral  motion  of  the  balloon  during  ascent, 
float,  or  descent.  This  is  likely  attributable  to  the  practical 
difficulty  of  incorporating  a  reliable  and  useful  wind  model. 
Also,  since  CRC  was  unable  to  locate  any  data  or  reference  to 
data,  it  is  reasonable  to  assume  that  only  in  rare  cases  is  a 
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balloon  subjected  to  low-speed  wind  tunnel  testing  to  formulate 
accurate  representations  of  a  balloon's  response  to  horizontal 
air  flow.  On  the  other  hand,  a  balloon's  vertical  motion  can  be 
understood  more  readily  since  the  vertical  component  of  motion  is 
more  repeatable  during  tests  and  actual  flights.  This  is  the 
apparent  explanation  for  the  presence  of  detailed,  validated 
models  of  a  balloon's  vertical  motion  and  for  the  absence  of  more 
comprehensive,  multi-dimensional  balloon  flight  performance 
models . 

The  search  for  an  atmospheric  model  for  the  simulation  of 
winds  in  BDPS  led  to  a  single  candidate,  NASA's  Global  Reference 
Atmosphere  Model  (GRAM).  GRAM  is  an  empirical  FORTRAN  computer 
simulation  of  the  earth's  atmosphere.  It  was  developed  by 
Georgia  Tech  under  contract  to  NASA's  Marshall  Space  Flight 
Center.  Dale  (ref.  6]  summarizes  the  capabilities  and  operation 
of  GRAM.  The  features  which  are  most  significant  to  BDPS 
development  are  that  GRAM  provides  a  worldwide,  12-month  database 
of  atmospheric  properties  including  wind  speed  and  wind  azimuth. 
The  atmospheric  properties  are  provided  in  statistical  form  (mean 
values  and  standard  deviations  from  the  mean)  in  order  to  best 
represent  the  actual  data  collected  over  several  decades.  NASA's 
most  significant  use  of  GRAM  was  in  the  aerothermal  design  of  the 
Space  Shuttle  Transportation  System  where  atmospheric  reentry 
heating  issues  had  to  be  resolved. 

Previously,  CRC  had  acquired  GRAM  and  had  become  familiar 
with  its  operation.  This  factor,  added  to  GRAM'S  suitability  to 
the  BDPS  mission,  made  it  an  obvious  choice  for  inclusion  in 
BDPS. 

2.2  MODEL  DEVELOPMENT  AND  INTEGRATION 

As  shown  in  the  Figure  1.0-1  schedule,  the  majority  of 
effort  for  this  contract  was  focused  on  modeling  taslcs.  First, 
atmosphere  models  and  balloon  models  were  implemented.  Then, 
those  models  were  integrated  with  a  vehicle  flight  dynamics 
program  to  develop  the  BDPS  capability. 
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2.2.1  Atmospheric  Model 

GRAM  is  a  combination  of  different  empirical  atmosphere 
models  corresponding  to  three  different  altitude  regimes.  Figure 
2. 2. 1-1  shows  these  regimes  and  the  models  used.  Between  0  and 
25  kilometers  atmospheric  pressure,  density  and  temperature  are 
simulated  by  the  4-Dimensional  World-Wide  Atmosphere  Model 
developed  by  Allied  Research  Associates  under  NASA  sponsorship. 
The  4-D  model  accesses  large  data  files  containing  empirically 
determined  atmospheric  parameter  profiles  for  a  grid  of  world¬ 
wide  locations.  A  modified  version  of  the  latitude-dependent 
model  of  G.  V.  Groves  is  used  for  altitudes  between  30  and  115 
kilometers.  This  model  was  modified  to  incorporate  longitudinal 
variations  in  pressure,  density,  and  temperature  in  addition  to 
the  latitudinal  variations.  For  altitudes  above  90  kilometers, 
the  model  of  L.  G.  Jacchia  is  used  to  compute  the  thermospheric- 
exospheric  values  of  atmospheric  pressure,  density,  and 
temperature.  Interpolation  between  the  4-D  and  the  Groves  values 
is  performed  between  25  and  30  kilometers,  and  a  fairing 
technique  is  used  between  the  Groves  and  the  Jacchia  models. 

Two  different  wind  models  are  used  in  GRAM.  The  geostrophic 
wind  model  is  applicable  for  altitudes  between  0  and  25 
kilometers  and  above  90  kilometers.  This  model  uses  finite 
differencing  of  pressure  values  to  estimate  horizontal  pressure 
gradients  which  are  then  used  to  compute  wind  velocities.  The 
wind  speed  model  for  the  middle  atmosphere  from  25  to  90 
kilometer  is  the  spherical  harmonic  wind  model  based  on  the 
spherical  harmonic  equation  which  is  a  function  of  latitude, 
longitude,  month,  and  altitude  and  whose  coefficients  are  fitted 
to  measured  data  and  vary  with  month  and  altitude.  Smooth 
fairing  between  spherical  harmonic  winds  and  geostrophic  winds  is 
performed  between  20  and  25  kilometers  and  between  90  and  95 
kilometers.  The  GRAM  wind  model  also  computes  values  for  mean 
vertical  winds  and  thermal  wind  shear.  The  mean  vertical  winds 
are  typically  of  the  order  of  a  centimeter/second  or  less  and 
represent  the  large-scale  mean  vertical  winds  characteristic  of 
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mean  meridional  circulation.  Wind  shear  values  are  calculated 
for  any  altitude  below  700  kilometers  but  are  reliable  only  below 
90  kilometers. 

Figure  2. 2. 1-2  shows  a  sample  wind  profile  provided  by  GRAM. 
The  figure  shows  the  east  wind  component  and  the  north  wind 
component  at  a  latitude  and  longitude  corresponding  to  Wallops 
Island,  Virginia.  The  January  database  was  used  by  GRAM  to 
produce  the  data  shown  in  this  figure.  The  graphs  show  the  mean 
wind  velocity  in  each  direction  along  with  the  mean  ±  3a  values. 
The  mean  ±  3o  curves  are  significant  because  99.7%  of  the  wind 
speed  values  will  fall  within  those  bounds. 


2.2.2  Balloon  Model 

CRC  decided  that  a  single,  simplified  model  could  be 
developed  to  model  both  zero-pressure  and  super-pressure  balloon 
aspects  which  influence  a  balloon's  drift  pattern.  The 
requirements  imposed  on  the  balloon  model  developed  were  that  it 
be  completely  configurable  through  input  data  rather  than  through 
coding  changes  and  that  it  provide  all  the  balloon-specific 
parameters  necessary  for  drift  pattern  production  to  the  vehicle 
flight  dynamics  portion  of  the  BDPS.  Figure  2. 2. 2-1  shows  a  top- 
level  flowchart  of  the  balloon  model,  and  the  actual  FORTRAN 
source  code  is  provided  as  Appendix  A.  As  shown  in  the  figure, 
the  only  modeling  difference  in  the  treatment  of  the  two  types  of 
balloons  is  that  a  zero-pressure  balloon  will  expel  excess  gas 
while  a  super-pressure  balloon  retains  the  excess  gas,  causing  an 
internal  pressure  greater  than  the  current  ambient  atmospheric 
pressure . 

By  making  the  balloon  model  data-conf igurable ,  the  source 
code  edit-compile-link  cycles  are  confined  to  code  development. 
This  feature  enhances  the  usefulness  of  the  balloon  model  for 
actual  production  of  drift  patterns  since  no  parameters  are 
"hard-wired"  into  the  source  code.  Table  2. 2. 2-1  lists  the  input 
values  required  to  configure  a  balloon  with  this  model.  The 
appropriate  units  for  the  input  values  are  also  shown  in  the 
table . 
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Table  2. 2. 2-1.  Balloon  Model  Input  Data 


Parameter 


Units 


Type  of  balloon 

(zero-  or  super-pressure) 

Gas  constant  of  lifting  gas 

Weight  of  lifting  gas 
Weight  of  payload  and  balloon 
Maximum  diameter  of  balloon 
Drag  coefficient 


meter^ 


second^ • “Kelvin 
newton 
newton 
meter 


The  model  imposes  spherical  geometry  on  the  balloons.  The 
"maximum  diameter"  input  parameter  in  Table  2. 2. 2-1  is  used  as  a 
spherical  diameter  when  determining  if  either  type  of  balloon  is 
fully  inflated.  Also,  the  aerodynamic  reference  area  for  fully 
or  partially  inflated  balloons  is  calculated  by  using  the 
assumption  that  the  balloon's  shape  is  always  approximately 
spherical.  Dwyer  (ref.  51  dwells  on  the  point  that  the  shape  of 
a  partially  inflated  balloon  must  be  treated  carefully  if  a 
simulation  is  to  attain  fidelity  with  the  actual  behavior  of  a 
balloon. 

The  balloon  model  provides  three  essential  output  quantities 
for  use  in  the  vehicle  flight  dynamics  program.  As  discussed 
above,  a  reference  area  is  calculated  and  serves  as  the  area  upon 
which  the  aerodynamic  drag  calculation  is  based.  Also,  the  model 
provides  the  magnitude  of  the  buoyant  force  which  is  equal  to  the 
weight  of  the  air  displaced  by  the  balloon.  Finally,  the  model 
reports  the  current  weight  of  lifting  gas  in  the  balloon  which 
can  only  change  if  a  zero-pressure  balloon  expels  excess  gas  as 
it  enters  its  float  altitude. 
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In  general,  CRC  recognizes  that  the  thermodynamic  fidelity 
of  the  balloon  model  described  in  this  report  suffers  from 
comparison  with  highly  detailed  models  such  as  those  of  Kreith 
(ref.  2]  or  Carlson  and  Horn  (ref.  3],  but  the  effort  required  to 
increase  the  fidelity  of  the  present  model  would  not  necessarily 
be  rewarded  by  significant  differences  in  the  resulting  simulated 
drift  patterns.  For  example,  the  addition  of  modeling  the 
thermal  stratification  of  the  lifting  gas  would  not  be  expected 
to  cause  a  noticeable  difference  in  the  distance  covered  by  a 
drifting  balloon  over  several  hours. 

2.2.3  Balloon  Flight  Dynamics/Kinematics  Model 

CRC's  proprietary  Multiple  Degree  of  Freedom  (MDOP) 
simulation  system  provided  the  trajectory  simulation  framework 
required  to  demonstrate  the  BDPS  capability.  MDOF  is  a  highly 
versatile  FORTRAN  computer  program  written  by  missile  design  and 
control  system  engineers  as  a  tool  for  analyzing  missile  flight 
and  engagement  performance.  It  is  designed  to  provide  rapid 
synthesis  of  a  wide  variety  of  air  vehicle  configurations  while 
allowing  direct  growth  of  simulation  fidelity  with  vehicle  design 
maturity.  MDOF's  modular  structure  and  program  control  features 
provide  the  user  with  the  essential  elements  of  all  flight 
simulations;  input/output,  equations  of  motion,  executive 
control/sequencing,  integration,  and  basic  configuration  modules. 

2.2.4  Balloon  Drift  Pattern  Simulation  (BDPS)  Integration 

The  integration  of  the  balloon  model  into  the  MDOF  framework 

was  simple  and  straightforward.  All  input  and  output  values  were 
passed  to  and  from  the  routine  through  the  respective  common 
blocks.  In  addition  to  the  three  main  outputs  discussed  above, 
the  balloon  model  also  made  available  calculations  of  free  lift 
percentage  and  internal  balloon  pressure.  Once  the  calculations 
were  made,  those  parameters  could  easily  be  added  or  deleted  from 
the  output  files  as  desired  by  the  user. 
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CRC  decided  to  execute  GRAM  off-line  and  reduce  its  output 
to  tabular  wind  data  which  could  then  be  inserted  into  an  MDOF 
data  file.  This  allows  for  a  single  wind  profile  (wind  speed  and 
wind  azimuth,  each  as  a  function  of  altitude)  to  be  generated  and 
then  used  for  the  entire  BDPS  execution.  The  wind  profile  for 
each  24-hour  flight  was  generated  for  the  balloon's  launch 
position.  The  use  of  a  single  wind  profile  throughout  a  24-hour 
flight  introduced  deviations  in  wind  speed  potentially  as  high  as 
15%  when  comparing  the  wind  speed  value  taken  from  the  input 
profile  to  the  value  which  GRAM  would  have  reported  for  the 
balloon's  current  position.  An  alternative  approach  would  have 
been  to  patch  the  entire  GRAM  model  directly  into  MDOF.  However, 
either  the  combined  program  would  be  terribly  inefficient  with 
GRAM  having  to  re-initialize  itself  during  each  MDOF  execution 
step,  or  GRAM  would  have  to  be  completely  restructured  and 
distributed  throughout  the  MDOF  code.  Neither  of  these  options 
were  feasible  under  the  SBIR  Phase  I  time  and  funding 
constraints.  A  sample  profile  as  used  by  BDPS  is  shown  in  Figure 
2. 2. 4-1,  where  the  wind  speed  values  have  been  calculated  by  a 
root-sum-square  combination  of  the  east  and  north  wind  components 
shown  previously  in  Figure  2. 2. 1-2. 

The  BDPS  data  flow  is  summarized  in  Figure  2. 2. 4-2.  Balloon 
design  data  is  combined  with  information  about  the  balloon's 
intended  mission.  That  data  is  used  to  build  an  input  data  file 
for  MDOF  and  provides  location  data  for  creating  a  wind  profile 
with  GRAM.  The  GRAM  wind  profile  is  then  inserted  into  the  MDOF 
input  data  file,  and  MDOF  is  executed  with  the  integrated  balloon 
model.  Depending  on  what  variables  were  requested  through  the 
input  file,  certain  variables  are  made  available  for  display 
using  CRC's  proprietary  PlotPlus  graphics  tool.  Among  other 
output  options,  PlotPlus  may  be  used  to  produce  balloon  drift 
patterns . 
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JANUARY:  WALLOPS  ISLAND,  VIRGINIA 


Wind  Profile 


BALLOON  DESIGN/ 
MISSION  DATA 


Figure  2.2A-2.  Data  Flow  for  Drift  Pattern  Production 
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2.2.5  Output  Format  Development 

After  establishing  contact  with  Lt.  Cdr.  Epley,  CRC 
discussed  DARPA's  needs  and  developed  sample  output  formats  which 
were  sent  to  Lt.  Cdr.  Epley  and  to  Les  Reading.  The  formats 
received  approval  and  are  used  in  this  report  to  present  the 
results  of  demonstrating  the  BDFS  capability. 

The  primary  output  format  requested  by  DARPA  was  a  plot 
showing  the  latitude  and  longitude  of  a  balloon  with  major  land 
masses  also  visible  for  reference.  Other  plots  of  interest  to 
DARPA  included  an  ascent  velocity  history  and  a  plot  of  the 
balloon's  altitude  versus  the  ground  range  from  the  balloon's 
launch  point.  Samples  of  these  two  output  formats  are  shown  in 
Figures  2. 2. 5-1  and  2. 2. 5-2. 

2.3  DRIFT  PATTERN  PRODUCTION  AND  ANALYSIS 

The  data  used  to  configure  the  zero-pressure  balloon  shown 
in  six  sample  flights  was  primarily  obtained  from  Les  Reading. 
Specifically,  he  provided  the  weight  of  the  balloon  fabric  and 
payload,  the  approximate  dimensions  of  the  fully  inflated 
balloon,  the  approximate  amount  of  lifting  gas  used  to  launch  the 
balloon,  and  the  type  of  lifting  gas.  He  also  provided  CRC  with 
estimates  of  the  initial  and  final  ascent  velocity  that  he  had 
observed  during  balloon  flight  tests. 

CRC  obtained  an  estimate  of  the  balloon  drag  coefficient 
from  Hoerner  (ref.  7].  The  coefficient  value  chosen  was  based  on 
the  balloon  description  given  CRC  by  Les  Reading,  who  described 
his  design  as  a  "pillow-shaped"  rectangular  balloon  with  the  long 
axis  oriented  vertically.  For  an  ascending  balloon,  such  a  shape 
could  be  most  closely  approximated  by  a  circular  cylinder  in 
axial  flow,  which  is  catalogued  by  Hoerner  as  having  a  pressure 
drag  coefficient  of  0.8.  This  seems  to  be  a  reasonable  choice 
since  the  use  of  this  number  produced  simulated  ascent  velocity 
profiles  which  resembled  those  reported  by  Les  Reading  from 
flight  tests.  It  should  be  noted,  however,  that  the  balloon 
model  computes  an  aerodynamic  reference  area  based  on  a  spherical 
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WALLOPS  ISLAND  LAUNCH 
HOURLY  POSITION  MARKS 


2ui.5-2.  Sample  Output:  Balloon  Altitude  vs.  Ground  Range 


volume,  for  which  Hoerner  reports  a  drag  coefficient  of  0.47. 
The  apparent  inconsistency  between  drag  coefficients  and 
reference  areas  serves  to  point  out  that  the  ascending  balloon's 
shape  and  drag  characteristics  are  not  modeled  to  a  high  level  of 
detail  in  this  approach.  In  spite  of  this,  the  model  produced 
results  which  track  well  with  the  reported  flight  test  results. 
The  simulated  ascent  rate  and  the  flight  test  ascent  rate  are 
both  shown  in  Figure  2.3-1. 

Lt.  Cdr.  Krai  identified  launch  sites  to  be  used  for  both 
zero-pressure  flights  and  for  a  super-pressure  flight.  The  six 
sites  from  which  zero-pressure  balloon  drift  patterns  are  to 
originate  are  listed  in  Table  2.3-1.  The  latitude  and  longitude 
values  for  each  location  were  estimated  by  CRC .  The  single 
location  to  be  used  for  a  super-pressure  launch  is  also  indicated 
in  the  table. 


For  each  of  the  six  launch  locations  listed  in  Table  2.3-1, 
four  plots  are  presented.  First,  an  overview  chart  shows  the 
balloon  drift  pattern  with  respect  to  nearby  land  masses. 
Second,  a  magnified  view  of  the  latitude-longitude  excursion  is 
shown.  The  third  plot  for  each  launch  shows  the  balloon's 
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SIMULATION 
FLIGHT  TEST  DATA 


Figure  2.3-1.  Comparison  of  Simulated  Ascent  Rate  and  Flight  Test  Ascent  Rate 


altitude  as  a  function  of  its  ground  range.  The  fourth  chart  for 
each  location  shows  the  balloon's  ascent  velocity  history  during 
the  time  of  ascent.  The  fifth  and  sixth  charts  for  each  launch 
site  show  the  January  speed  and  azimuth  wind  profiles  which  were 
used  to  generate  that  drift  pattern.  This  series  of  plots  can  be 
found  in  Figures  2.3-2  through  2.3-37. 

Figure  2.3-38  shows  several  24-hour  drift  patterns  for 
balloons  launched  from  Wallops  Island.  For  this  figure,  the  BDPS 
was  provided  with  several  wind  profiles  which  represent  the 
statistical  variations  to  be  expected  during  January.  Bounded  by 
the  drift  patterns  using  the  mean  +  3o  winds,  the  area  over  which 
the  balloons  drifted  represents  99.7%  of  the  area  over  which  the 
balloons  could  conceivably  drift  for  that  24-hour  period  in 
January . 

For  the  Wallops  Island  launch  site,  four  drift  patterns  were 
produced  showing  the  effect  of  winds  during  different  months  of 
the  year.  Figure  2,3-39  shows  the  four  drift  patterns,  produced 
using  January,  April,  July,  and  October  winds. 

Figure  2.3-40  through  2.3-67  show  the  drift  pattern  of  the 
super-pressure  balloon  launched  from  Adak  Island  and  allowed  to 
drift  at  120,000  feet  for  one  year.  CRC  was  given  latitude  in 
the  configuration  specification  of  this  balloon  since  the  only 
available  information  was  its  design  altitude  and  its  50-lb 
payload.  The  drift  pattern  was  segmented  among  several  pages  for 
clarity  of  presentation. 
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.  Altitude  Versus  Ground  Range  for  24-Hour  Wallops  island  Flight 
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24-Hour  Drift  Pattern  for  Vandenberg  AFB  Flight 


24>Hour  Drift  Pattern  for  Vandenberg  AFB  Flight 
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Altitude  Versus  Ground  Range  for  24-Hour  Vandenberg  AFB  Flight 


Ascent  Velocity  History  for  24-Hour  Vandenberg  AFB  Flight 
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Figure  2.^12.  Wind  Speed  Profile  at  Vandenberg  AFB  in  January 
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Figure  2.3-16.  Altitude  Versus  Ground  Range  for  24-Hour  Adak  Island  Flight 
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Figure  2.3-20.  24-Hour  Drift  Pattern  for  Fiight  North  of  Puerto  Rico 
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Figure  2.3-21 .  24-Hour  Drift 
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Figure  2.3-22.  Altitude  Versus  Ground  Range  for  24-Hour  Flight  North  of  Puerto  Rico 
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Figure  2.3-28.  Altitude  Versus  Ground  Range  for  24-Hour  Mediterranean  Sea  Flight 


Ascent  Velocity  History  for  24>Hour  Mediterranean  Sea  Flight 
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Figure  2.3-33.  24-Hour  Drift  Pattern  for  indian  Ocean  Fiight 


Altitude  Versus  Ground  Range  for  24-Hour  Indian  Ocean  Flight 
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Figure  2.3-36.  Wind  Speed  Profile  in 
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Figure  2.3-37.  Wind  Azimuth  Profile  in  Indian  Ocean  in  January 


lur  Drift  Patterns  for  Wallops  Island  Flight  with 
Statistical  January  Winds 


Hour  Drift  Patterns  for  Wallops  Island  Flights  Showing 
Effects  of  Seasonal  Wind  Variations 
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Figure  2.3-40.  Super-Pressure  Balloon  Drift  Pattern  at  120,000  ft  Altitude: 
Day  1  through  Day  120  (Jan  1  —  Apr.  30) 
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Pressure  Balloon  Drift  Pattern  at  120,000  ft  Altitude: 
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Figure  2.3-42.  Super-Pressure  Balloon  Drift  Pattern  at  120,000  ft  Altitude: 
Day  120  through  Day  158  (Apr.  30  —  Jun.  7) 
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Figure  2.3-43.  Super-Pressure  Balloon  Drift  Pattern  at  120,000  ft  Altitude: 
Day  120  through  Day  158  (Apr.  30  —  Jun.  7) 
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Figure  2.S-44.  Super-Pressure  Balloon  Drift  Pattern  at  120,000  ft  Altitude: 
Day  158  through  Day  173  (Jun.  7  —  Jun.  22) 


e  2.3-45.  Super-Pressure  Balloon  Drift  Pattern  at  120,000  ft  Altitude: 
Day  158  through  Day  173  (Jun.  7  —  Jun.  22) 
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Figure  2.3-46.  Super-Pressure  Balloon  Drift  Pattern  at  120,000  ft  Altitude: 
Day  173  through  Day  187  (Jun.  22—  Jul.  6) 


68 


Figure  2.^-47.  Super-Pressure  Balloon  Drift  Pattern  at  120,000  ft  Altitude 
Day  173  through  Day  187  (Jun.  22—*  Jul.  6) 
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e  2.3-48.  Super-Pressure  Balloon  Drift  Pattern  at  120,000  ft  Altitude: 
Day  187  through  Day  203  (Jul.  6—  Jul.  22) 
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Super-Pressure  Bafioon  Drift  Pattern  at  120,000  ft  Altitude: 
Day  203  through  Day  223  (Jui.  22~  Aug  11) 
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3uper-Pressure  Balloon  Drift  Pattern  at  120,000  ft  Altitude: 
Day  203  through  Day  223  (Jul.  22—  Aug  11) 
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Pressure  Balloon  Drift  Pattern  at  120,000  ft  Altitude: 


Super-Pressure  Balloon  Drift  Pattern  at  120,000  ft  Altitude: 
Day  232  through  Day  237  (Aug  20  —  Aug  25) 


9r-Pressure  Balloon  Drift  Pattern  at  120,000  ft  Altitude: 
232  through  Day  237  (Aug  20  —  Aug  25) 


Pressure  Balloon  Drift  Pattern  at  120,000  ft  Altitude: 


; _ I _ I _ I _ I — I — I _ I _ I _ I _ I _ I _ L 

o  o  o 

ID  ID 

I 


930  BOniliVI 


o 

o 

OJ 


o 

ID 


O 

CO 


O 
O  O 
O 

^  1 


I 


0) 

■o 

3 


< 

ti. 

o 

o 

o 

o' 

CM 

T” 

(0 

c 

k. 

o 

s 

a 


c 

o 

o 

0 

CD 

£ 

3 

<A 

M 

O 

k> 

a. 

ii. 

o 

a 

3 

(0 


In 


CM 

A 

h. 

3 

CD 

u! 


78 


Day  237  through  Day  253  (Aug  25  ~  Sep  10) 


253  _  _  263  26B 
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Figure  2.3-58.  Super-Pressure  Balloon  Drift  Pattern  at  120,000  ft  Altitude: 
Day  253  through  Day  268  (Sep  10  —  Sep  25) 
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Pressure  Balloon  Drift  Pattern  at  120,000  ft  Altitude: 
>3  through  Day  268  (Sep  10  —  Sep  25) 


930  3aniiiv~i 


Day  268  through  Day  283  (Sep  25  —  Oct  10) 


Pressure  Balloon  Drift  Pattern  at  120,000  ft  Altitude 
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Figure  2.3-62.  Super*Pres8ure  Balloon  Drift  Pattern  at  120,000  ft  Altitude: 
Day  283  through  Day  310  (Oct  10  —  Nov  6) 
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Pressure  Balloon  Drift  Pattern  at  120,000  ft  Altitude: 


Figure  2.3-64.  Super-Pressure  Balloon  Drift  Pattern  at  120,000  ft  Altitude: 
Day  310  through  Day  332  (Nov  6  —  Nov  28) 
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Super-Pressure  Balloon  Drift  Pattern  at  120,000  ft  Altitude: 
Day  332  through  Day  365  (Nov  28  —  Dec  31) 
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Pressure  Balloon  Drift  Pattern  at  120,000  ft  Altitude: 


3.0  CONCLUSIONS 

CRC  has  demonstrated  a  capability  to  produce  simulated 
balloon  drift  patterns  using  digital  computer  simulation 
techniques.  Where  possible,  available  balloon  flight  test  data 
was  used  for  comparison  with  the  simulation  results,  showing 
similar  trends  with  reasonable  agreement. 


The  models  used  and  developed  for  this  effort  proved 
adequate  but  not  optimal.  A  particular  need  exists  for  a  closer 
coupling  between  GRAM  and  MDOF  in  order  to  reduce  potential 
inaccuracies  which  may  arise  from  using  a  single  wind  profile 
over  the  extent  of  a  balloon's  drift.  More  tightly  integrated 
BPDS  models  would  also  serve  to  reduce  the  effort  associated  with 
the  production  of  a  single  drift  pattern,  whether  it  is  a  24-hour 
or  1-year  mission. 
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APPENDIX  A 


BALLOON  MODEX  EXXm»N  SOURCE  CODE 


A-1 


SUBROUTINE  BUOY 

002  C 
004  C 

005  C  REVISIC»I  DATE  -  16-JAN-1989 

006  C  - 

007  C  PURPOSE 

008  C  ,  BUOYANT  FXDRCE  MODEL 

009  C 

010  C  PROGRAM  ROUTINES 

Oil  C  .  NCttJE 

012  C 

013  C  UTILITY  ROUTINES 

014  C  .  TLUlEI  -  SINGLE  VARIABLE  TABLE  LOOK-UP 

015  C  .  VMAG3  -  VECTOR  (3)  MACaJITUDE 

016  C 

017  C  SYSTEM  LIBRARIES 

018  C  .  (NCWE) 

019  C 

021  C 

022  DIMENSIOJ  ID  (  9999  ) 

023  C 

024  PARAMETER  (  AIRMW  -  28.97  ) 

025  C 

026  C: ::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::  :C 
027  C 

028  C**  SIMULATION  CC»lPUTATIONAL  COMMON  *  *  *  * 

029  C 

030  COMMCXa  /  DCCW  /  D  (  9999  ) ,  V  (  5000  ),  C  (  100  ) ,  JD  (  800  ) 

031  C 

032  C**  UPDATE  OOMMC»4  ************** 

033  C 

034  COMMCW  /  UCOM  /  UP  (  200  ) 

035  C 
036  C 
037  C 
038 
039  C 
040 
041 
042  C 
043 
044 
045 
046 
047 
048 
049 
050 
051 
052 


EQUIVALENCE 

( 

D( 

EQUIVALENCE 

( 

D( 

EQUIVALENCE 

( 

D( 

EQUIVALENCE 

( 

V( 

EQUIVALENCE 

( 

V( 

EQUIVALENCE 

( 

V( 

EQUIVALENCE 

( 

V( 

EQUIVALENCE 

( 

V( 

EQUIVALENCE 

( 

V( 

EQUIVALENCE 

( 

V( 

EQUIVALENCE 

( 

V( 

EQUIVALENCE 

( 

V( 

EQUIVALENCE 

( 

V( 

1) ,ID(1)  ) 

66),RGAS  ),( 
69),GSFACT) 

2) ,DT  ),( 

51),GXI  ),( 

69), PRESS  ),( 
96),AREF  ),( 

126), QS  ),( 

243),XDWT  ),( 
246),FBUDYX),( 
245),FBUOYT) 
267),BPRESS),( 
270),FLIFT  ) 


D(  67),DIAMAX),( 


V(  5), TABS  ),( 
V(  52),GYI  ),( 

V(  70), DENS  ),( 
V(  101),TnATA  ),( 
V(  214),TSUBA  ) 

V(  244),XDlWr  ) 

V(  247),FBUOYY),( 

V(  268),SPDIAM),( 


D(  68),GASnN  ) 


V(  7),X3600  ) 
V(  53),GZI  ) 

V(  71),DYNP  ) 
V(  115),TTLU  ) 


V(  248),FBUOYZ) 
V(  269),GASWr  ) 


nononooooooononn 


c 


053 


054 

EQUIVALENCE 

(  C( 

1),PI  ),( 

C( 

65), GO 

) 

055 

EQUIVALENCE 

(  C( 

81),GMASS  ),( 

C( 

82),DEXPEL),( 

C( 

83), VOLUME) 

056 

EQUIVALENCE 

(  C( 

84 ), VOLMAX) ,  ( 

C( 

85), GIMAG 

),( 

C( 

86),TMP6  ) 

057 

EQUIVALENCE 

(  C( 

87),  WEIGHT),  ( 

C( 

88), WGHTO 

),( 

C( 

89),GASOUT) 

058 

059  C 

EQUIVALENCE 

(  C( 

90),GMAX  ) 

060 

061  C 

EQUIVALENCE 

(ID( 

2), IBUOY  ) 

062 

EQUIVALENCE 

(JD( 

55),IRWGT  ),(JD( 

438) , ICWGT 

) 

063  C 

064  SAVE 

065  C 

066  C - C 

067  C 

068  C  CALCULATE  MAGNITUDE  OF  INERTIAL  GRAVITY  VECTOR 
069  C 

070  CALL  VMAG3  (  GXI,  GIMAG,  TMP6  ) 

071  C 

072  C  FIND  VOLUME  OF  LIFTING  G?^  AT  CURRENT  AMBIENT  PRESSURE  AND 
073  C  TEMPERATURE  USING  IDEAL  GAS  LAW:  V  =  MRT/P 
074  C 

075  OlASS  =  GASWT  /  GO 

076  VOLUME  -  (21ASS  *  RGAS  *  TSUBA  /  PRESS 

077  C 

078  C  SET  INTERNAL  BALLOON  PRESSURE  EQUAL  TO  AMBIENT  PRESSURE 
079  C 

080  BPRESS  =  PRESS 

081  C 

082  C  CALCULATE  CURRENT  DIAMETER  ASSUMING  SPHERICAL  SHAPE 
083  C 

084  SPDIAM  -  (  6.0  *  VOLUME  /  PI  )  **  0.333333 

085  C 

086  C  IBUOY*!  ->  ZERO-PRESSURE  BALLOCW 

087  C 

088  IF  ( IBUOY. BQ.l)  THEN 

089  C 

090  C  CALCULATE  FREE  LIFT  FOR  ZERO-PRESSURE  BALLOON;  FIRST, 

091  C  COMPUTE  INDEPENDENT  VARIABLE  FOR  TABLE  LOOK  UP  TO  GET 

092  C  CURRENT  WEIGHT  EXCLUDING  WEIOTT  OF  LIFTING  GAS 

093  C 

094  TTLU  =  TABS  -  TDATA 

095  TMP6  =  TTLU*XDIWT 

096  CALL  TLUlEI  (TMP6,  D(IRWC?r),  ICWGT,  WGHTO) 

097  WEIGHT  -  W3HTO*XDWT/GO 

098  C 

099  FLIFT  -  100.0  *  ((OlASS  *  (AIRMW/GASMW  -  1.0))  /  WEIGHT  -1.0  ) 

100  C 

101  C  COMPARE  TOE  IDEAL-GAS/SPHERICAL  DIAMETER  TO  THE  MAXIMUM 

102  C  EFFECTIVE  SPHERICAL  DUtfETER 

103  C 

104  IF  ( SPDIAM. GT.DIAMAX)  THEN 

105  VOLMAX  -  (  PI  /  6.0  )  *  DIAMAX**3 

106  aiAX  -  (PRESS  *  VOLMAX)  /  (RGAS  *  TSUBA) 


IF  (GMASS.GT.(GMAX*GSFACT) )  THEN 


L07  C 
LOS 
L09  C 

LIO  C  CALCULATE  THE  MASS  EXPULSION  RATE  FOR  EXCESS  LIFTING  GAS 

Lll  C 

L12  DEXPEL  =  (GMASS  -  GMAX*GSFACT)  /  (DT  *  8.0) 

L13  GASOUT  =  DEXPEL  *  DT 

L14  GMASS  =  GMASS  -  GASOUT 

L15  GASWT  =  GMASS  *  GO 

L16  ENDIF 

L17  C 

L18  VOLUME  -  VOLMAX 

L19  BPRESS  -  GMASS  *  RGAS  *  TSUBA  /  VOLUME 

L20  EJ4DIF 

L21  C 

L22  C  IBUOY«2  ->  SUPERPRESSURE  BALLOON 

L23  C 

L24  ELSE  IF  (IBUOY.EQ.2)  THEN 

L25  C 

L26  C  COMPARE  THE  IDEAL-GAS/SPHERICAL  DIAMETER  TO  THE  MAXIMUM 

L27  C  EFFECTIVE  SPHERICAL  DIAMETER 

L28  C 

L29  IF  (SPDIAM.GT.DIAMAX)  THEN 

L30  SPDIAM  -  DIAMAX 

L31  VOLUME  «=  (  PI  /  6.0  )  *  SPDIAM**3 

L32  BPRESS  =  GMASS  *  RGAS  *  TSUBA  /  VOLUME 

L33  ENDIF 

L34  C 

L35  ENDIF 

L36  C 

L37  C  COMPUTE  TOTAL  BUOYANT  FORCE  (WEIGHT  OF  THE  DISPLACED  FLUID) 

L38  C 

L39  FBUOYT  -  VOLUME  *  DENS  *  GIMAG 

L40  C 

L41  C  COMPUTE  AERCX)YNAMIC  REFERENCE  AREA  AND  QS  PRODUCT 

L42  C 

L43  AREF  -  (  PI  /  4.0  )  *  SPDIAM**2 

L44  QS  -  DYNP  *  AREF 

L45  C 

L46  C  DECOMPOSE  TOTAL  BUOYANT  FORCE  INTO  BODY-REFEREJJCE  COMPC»IENTS  GIVEN 
L47  C  THAT  BALLOW  X-AXIS  IS  FORCED  TO  BE  ORIENTED  ALONG  GRAVITY  VECTOR 
L48  C 

L49  FBUOYX  -  FBUOYT 

L50  FBUOYY  -  0.0 

L51  FBUOYZ  »  0.0 

L52  C 

L53  RETURN 

L54  END 


